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Abstract. - A continuum version of the vortex-boson duality in (3+1) dimensions is formulated 
and its implications studied in the context of a pair Wigner crystal in underdoped cuprate super- 
conductors. The dual theory to a phase fluctuating superconductor (or superfluid) is shown to 
be a theory of bosonic strings interacting through a Kalb-Ramond rank-2 tensorial gauge field. 
String condensation produces Higgs mass for the gauge field and the expected Wigner crystal 
emerges as an interesting space-time analog of the Abrikosov lattice. 
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When particles in quantum many body systems inter- 
act very strongly standard perturbation techniques break 
down and, in dimensions greater than 1, dualities often 
provide the only insights into the physics of such systems. 
Duality transformations, in general, map the strongly cou- 
pled sector of one theory onto the weakly coupled sector 
of another. The original Kramers- Wanier duality [1] for 
the Ising ferromagnet represents a prime example of such 
a mapping. Dualities permeate modern statistical, con- 
densed matter and particle physics, and have emerged re- 
cently as a key tool in string theory. 

In condensed matter physics perhaps the most useful 
and influential duality is the one connecting vortices and 
bosons in two spatial dimensions [2-4] . This duality, here- 
after referred to as Lee-Fisher duality, maps the system 
of interacting bosons in (2+l)D onto a fictitious super- 
conductor in an external magnetic field whose flux in the 
temporal direction is proportional to the density of the 
original bosons. It shows that Mott insulator, proximate 
to the phase fluctuating boson condensate, can be viewed 
as the Abrikosov vortex lattice of the dual superconductor. 
This deep connection has been exploited in modeling sys- 
tems ranging from quantum spins to fractional quantum 
Hall effect, and most recently cuprate superconductors. 

In cuprates such considerations are motivated by the ex- 
perimental findings of static checkerboard patterns in the 
charge density of very underdoped samples [5-7] which 
have been interpreted as evidence for a Cooper pair 
Wigner crystal (PWC) [8]. The latter can be most nat- 
urally understood by appealing to the Lee-Fisher duality 
[9-11]. However, recent analysis of the vibrational modes 



of such a PWC [12, 13] indicates that it is 3-dimensional 
(in the sense that vibrations propagate in all 3 space di- 
mensions) and it is thus unclear how the inherently two- 
dimensional Lee-Fisher duality applies to this situation. 
The problem can be stated as follows. A key role in the 
formulation of the Lee-Fisher duality is played by vor- 
tices which appear (in pairs of opposite vorticity) near 
the transition to the Mott insulating phase as quantum 
fluctuations of the system. The dual relationship between 
bosons and vortices however exists only in two spatial di- 
mensions where the latter can be regarded as point par- 
ticles. In three space dimensions vortices form oriented 
loops and can no longer be thought of as particles. The 
question thus arises how to understand the formation of 
a PWC in the three dimensional phase fluctuating super- 
conductor indicated by experiments [5-7]. 

In this Letter we point a way out of this conundrum 
by constructing a (3+1) dimensional implementation of 
vortex-boson duality using a representation of vortex loops 
as relativistic bosonic strings. We then show that string 
condensation indeed produces a ground state that can be 
characterized as an insulating crystal of Cooper pairs and 
discuss some of its unique properties. 

We remark that the lattice formulation of such (3+l)D 
duality has been given long time ago [14, 15] and was used 
recently to study exotic fractionalized phases [16, 17] in 
(3+l)D. Here, by contrast, we formulate a continuum 
version which shows that a boson (Cooper pair) crystal 
can emerge from a phase fluctuating superfluid (super- 
conductor) even in the absence of any underlying lattice 
structure. This is exactly the limit apparently relevant 
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to cuprates [12]. On the formal side this continuum ap- 
proach also reveals an intimate connection to the theory 
of bosonic strings and enables us to employ in our calcu- 
lations some of the string theory technology. 

We begin by considering a continuum theory of 
a superconductor in (3+1) space-time dimensions de- 
fined by the Euclidean partition function Z = 
J £>[*, **] exp(- /(f dr J d 3 x£) with * = |*|e je a scalar 
order parameter and the Lagrangian density 



c = \k\{d. 



2ieA ll )*\ 2 + Um 2 ). 



(1) 



The Greek index fj, — 0,1,2,3 labels the temporal and 
spatial components of (3+1) dimensional vectors, and we 
use natural units with Ti = c = 1. U is a, potential func- 
tion that sets the value of the order parameter \I/ in the 
superconducting state in the absence of fluctuations. The 
electromagnetic vector potential A is explicitly displayed 
in order to track the charge content of various fields. If we 
allowed A to fluctuate then Eq. (1) would coincide with 
the well known Abelian Higgs model. 

We now focus on the fluctuations in the phase 8 by fixing 
the amplitude \^\ = "fo at the minimum of U, 



(2) 



where K = K^q represents the phase stiffness. The 
first few steps of the duality mapping proceed just as 
in (2+l)D. We first decouple the quadratic term with 
a real auxiliary field, W^, using the familiar Hubbard- 
Stratonovich transformation, obtaining 

C = JK W ^ + iW ^ Q - 2eA ^ + iW »(W- ( 3 ) 

We have also decomposed the phase into a smooth part 9 S 
and singular part 9 containing vortex lines. 

Gaussian integration over 6 S leads to a constraint 



0. 



(4) 



which reflects conservation of electric charge. In (2+l)D 
one enforces this constraint by expressing as a curl of 
an auxiliary gauge field. The curl operation, however, is 
meaningful only in 3 dimensions and herein lies the diffi- 
culty with higher dimensional duality. In (3+l)D we may 
enforce the constraint (4) by writing 



W li = e 



fj,i/ctj3 



d v B, 



a/3, 
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where e^ap is the totally antisymmetric tensor and B a $ is 
antisymmetric rank-2 tensor gauge field. Substituting Eq. 
(5) back into the Lagrangian and performing integration 
by parts in the term containing O we obtain 



3K 



(6) 



d a B{jy + dpB ia + d^B a fj is the tensorial 



where £f Q/ 3 7 

field strength which should be thought of as a generaliza- 
tion of the Maxwell field strength F^ v . 

The above Lagrangian exhibits several notable features. 
First, it possesses invariance under the gauge transforma- 
tion 
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for an arbitrary smooth vector function A M . The square 
brackets represent antisymmetrization, e.g. d[ a A^ = 
d a hp — dph a . This gauge invariance reflects conserva- 
tion of vorticity in the original model. Second, from the 
last term in C we may immediately deduce that the elec- 
tric four-current is related to B by j M = 2e(e fJil y a /}d u B a p). 
The charge density, in particular, can be written as 



p = jo = 2e(e ijk diB jk ), 



(8) 



where Roman indices run over spatial components only. 

The second term in C informs us that field B is mini- 
mally coupled to the "vortex loop current" 



(9) 



which is an antisymmetric rank-2 tensor quantity. For 
a smooth function the right hand side of Eq. (9) would 
vanish since the derivatives would commute. In 3 spatial 
dimensions, single valuedness of e 1 ®^ permits line singu- 
larities in 0(x) such that it varies by an integer multiple 
of 2ir along any line that encircles the singularity. These 
are the vortex loops. 

To implement the duality transformation we now shift 
our point of view from the phase field <d(x) to the vor- 
tex loop worldsheets. These describe the evolution of vor- 
tex loops in imaginary time and should be thought of in 
analogy with worldlines of point particles. The world- 
sheets are specified by a set of 2-parameter vector func- 
tions (0-1,0-2) where n labels the individual loops. We 
take G\ to be time-like, and correspondingly vary between 
and the inverse temperature (3, and 02 space-like, which 
by convention varies from to 2tt for closed loops. (Since 
vortices can only terminate on magnetic monopoles we 
consider only closed vortex loops here.) Clearly, given a 

set of worldsheets one can reconstruct the phase field 
0(x) up to any smooth contribution. 

A surface element of a worldsheet is characterized by a 
rank-2 antisymmetric tensor 
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It is straightforward to show that the loop current (9) is 
related to the worldsheet by 



a^(x) = 2nJ2 [ d 2 <J^5 (x™ - x) . (11 



This relation allows us to rewrite the partition function 
as a functional integral over the vortex loop worldsheets 
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xjTl We thus have Z = JV{X] exp(-S) with 



d xH nR ^ + 5i nt + 5j ac . 



(12) 



We recognize the first line as the celebrated Nambu-Goto 
action [18] for bosonic strings propagating in the pres- 
ence of a background Kalb-Ramond gauge field B^ v [19]. 
The first term can be interpreted as the reparametriza- 
tion invariant surface area of the string worldsheet with 
the string tension T. Although such term does not ap- 
pear explicitly in Eq. (6) it would arise in a more careful 
treatment of our starting theory (1) had we retained the 
cost of the suppression of the order parameter amplitude 
1^1 near the vortex core. The second and the third terms 
follow directly from Eq. (6) and describe long range inter- 
actions between strings mediated by the superflow, now 
represented by the Kalb-Ramond gauge field. 

<Sint contains short range interactions between strings 
that would also arise from a more careful treatment of 
the core physics. Finally 5j ac represents the Jacobian of 
the transformation from phase variable to string world- 
sheets X^\ This last term plays an important role in 
the quantization of our string theory. It is well known 
that a fundamental string can be consistently quantized 
only in the critical dimension, which for bosonic string is 
D = 26 [20] . A question then arises as to how we quantize 
our vortex strings in (3+l)D; after all we started from a 
well defined field theory (1) and we expect the string the- 
ory (12) derived from it to also be well behaved. The an- 
swer lies in the fact that our strings are not fundamental; 
rather they are Nielsen-Olesen type strings [21] with in- 
trinsic thickness defined by the core size. It was shown by 
Polchinski and Stromingcr [22] that terms in <Sj ac , which 
would be absent in the case of a fundamental string, pre- 
cisely cancel the conformal anomaly responsible for the 
high critical dimension. Vortex strings are indeed well be- 
haved in the physical dimension. 

We are now ready to complete the duality mapping. 
Our main goal will be to understand the string-condensed 
phase, analogous to the vortex-condensed phase in the 
(2+l)D Lee-Fisher duality. To this end we must pass 
to second quantized string theory, a "string-field the- 
ory" . This can be done rigorously in the so called light- 
cone gauge [23] or using the Becchi-Rouet-Stora-Tyutin 
(BRST) procedure [24]. Here we opt for a less rigorous 
but physically much more transparent procedure devised 
in Ref. [25] which provides a straightforward route towards 
the description of the string-condensed phase. 

The central concept in the string-field description is 
the wave functional $[X], a complex- valued functional 
defined on the space of one-parameter string trajectories 
{X^fa), 02 = (0, 27r)}, which we regard as cross sections 
of the string worldsheet X fl (ai,a2) at fixed value of o\. 
The physical significance of &[X] is most readily visualized 



in the so called static parametrization: in a chosen Lorentz 
frame of reference take X = a\ = t and X = X(r, a), 
with t the imaginary time. <&[t, X] then represents the 
quantum mechanical amplitude for finding the string in 
configuration {X(ct),ct = (0, 2ir)} at time r. 

The second quantized action for the string functional 
takes the form [25] 

S = J V[X] J daVh[\{5/5Z^-2TnB^)<5>[X]\ 2 

+ Ml s \^[X]\ 2 } + ^ j 'd 4 xH 2 a/37 +Si nt . (13) 

The plaquette derivative S/dT,^ quantifies the variation 
of the functional &[X] upon modifying the path X by an 
infinitesimal loop AX which sweeps the surface element 
ST,^. The loop space metric h — (dX^{a) / da) 2 is needed 
to preserve the reparametrization invariance of the action. 
S[ nt represents short range string interactions and contains 
terms cubic and higher order in |$|. All contributions 
quadratic in |$| have been folded into the effective string 
mass A^ ff. The action (13) remains invariant under the 
gauge transformation (7) if we require $ to transform as 



$[A] <f>[A]e 



-4-rri J dX^A^ 



(14) 



String condensation occurs when M 2 S becomes nega- 
tive. The string functional then develops nonzero vacuum 
expectation value, (0|$[A]|0) ^ 0. The simplest case is 
that of a uniform string condensate, 



(0|$[A]|0) = $ = const. 



(15) 



Physically, this simply means that any string configura- 
tion is equally probable. This ansatz, however, cannot 
describe a phase disordered superconductor. To see this 
note that substituting Eq. (15) into action (13) produces 
a mass term for the Kalb-Ramond gauge field. Such a 
mass term then leads to the Mcissner effect: the gauge 
field is expelled from the interior of the sample, = 0. 
In view of Eq. (8), this corresponds to complete expulsion 
of charge from the system, which is not the situation we 
are interested in. 

What we seek is the analog of the Abrikosov vortex state 
in which the field can penetrate in quantized increments. 
We thus consider a more general ansatz, which allows both 
the amplitude and the phase of &[X] to vary: 

(0|$[X]|0) = $ oe /^[CV^ln/ W +2 7 r^.n,(X)]_ 

Here X = X(a), X' = d a X(a), f{x) and Q^x) are 
real scalar and vector functions parametrizing the func- 
tional, and C is a parameter with the dimension of inverse 
length. f(x) is nonnegative and should be thought of as 
the space-time-dependent amplitude of the string conden- 
sate. Specifically, / = 1 corresponds to uniform conden- 
sate amplitude while / > 1 (/ < 1) describes its local 
enhancement (depletion), ^(x) determines the phase of 
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the string condensate. Substituting Eq. (16) to (13) we 
obtain S = J d 4 x£ with 



2 
1 



C = -^[K 2 .f 2 (d [ ,n ij] -2B, v ) 2 +( 2 (djr + V(f 2 )] 

(17) 



+ 3K H <*Pr 

We observe that any smooth part of f2 can be eliminated 
by the gauge transformation (14). Thus, only the singular 
part of £1 has physical significance. Indeed we note that 
it is permissible for £1 to be multiply valued as long as the 
wave functional (16) remains single valued. A configura- 
tion of specific interest to us contains mono-poles in the 
spatial part of Cl = (Sl n , fi), 



V.(Vxft)=^Q/)(x-4 



(18) 



where x a and Q a label the position and the charge of the 
a-th monopolc. Single valucdncss of (16) demands that Q a 
be integer. We shall see that such singularities represent 
sources for B^, just as vortices in a superconductor act 
as sources for the magnetic field. 

We now analyze the action (17) in the presence of static 
monopole configurations in tt. To this end we adopt a 
dual mean-field approximation (DMFA) which neglects 
quantum fluctuations of all the fields. We emphasize 
that in terms of the original phase degrees of freedom, 
DMFA describes a highly nontrivial quantum fluctuating 
state. In addition, we perform a dual "London" approx- 
imation, /(x) = 1, which should be adequate as long as 
the monopoles are relatively dilute. (This approximation 
fails in the small region near the monopole center where 
/ — ► 0.) The ground state energy of the system can then 
be written as 



7r <& (9[jfl,] - 2B lj ) 2 + —(cijkdiBjkY 



(19) 

Minimizing with respect to Bij leads to the Euler- 
Lagrange equation 



1 



2eK 



tijkdkP — 0, 



(20) 



where we used Eq. (8). Next, acting on all terms by e%jid\ 
and defining a dual "penetration depth" A^ 2 = 2tt 2 ^qK, 
we obtain an equation for charge density p(x) 



p-X 2 d \7 2 p = 2eV • (V x SI). 



(21) 



This equation resembles the London equation for the 
^-component of magnetic field in the presence of an 
Abrikosov lattice of vortices and can be analyzed by simi- 
lar methods. The key difference is that, in light of Eq. (18), 
the right hand side describes a collection of point sources 
in three space dimensions whereas Abrikosov vortices are 
line singularities described by <5' 2 ). Below we briefly sum- 
marize some main results of this analysis and the relevant 
details will be given elsewhere [26]. 



Eq. (21) can be solved for an arbitrary arrangement of 
monopole positions and charges to obtain 



p(x) = 2eJ2Qa 



g-|x-x a |/A<j 

47rA 2 |x-x a | 



(22) 



It is easy to show that the total electric charge associ- 
ated with a monopole is 2eQ a ; the charge is quantized 
in the units of 2e, as expected. At finite charge density, 
monopoles with like charges repel by Yukawa potential 
<~ er r l Xd jr and the ground state is a Bravais lattice of 
elementary (Q a = 1) monopoles. This leads to periodic 
modulation in p(x) with charge 2e per unit cell: a pair 
Wigner crystal in three space dimensions. 

An appealing overall picture thus emerges. Vortex 
loops in a (3+l)-dimensional superconductor (or super- 
fluid) can be efficiently described as bosonic strings in- 
teracting through a rank-2 tensorial Kalb-Ramond gauge 
field B^ v . In the non-superconducting phase, strings pro- 
liferate and condense, producing Higgs mass for the gauge 
field. In the Higgs phase the only way for B^ v to pene- 
trate into the bulk of the system is to set up quantized 
monopolc-likc singularities in the phase Q.^ of the string 
condensate wave functional. These singularities then act 
as point sources for B^ v . The associated electric charge 
density p, which is closely related to the Kalb-Ramond 
field strength H^x, is then governed by a London- like 
equation (21). For a periodic array of point sources, such 
as will form at finite charge density, a 3-dimensional pair 
Wigner crystal emerges with charge distribution given by 
Eq. (22). 

The duality discussed above establishes vortex-loop con- 
densation as a concrete mechanism for the formation of 
a pair Wigner crystal in a 3-dimensional quantum phase 
fluctuating superconductor. It explains how a 3d PWC 
can form in underdoped cuprates and allows for detailed 
computations of its structure and vibrational modes [26] 
which are of direct experimental interest. 
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